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NON-MAXWELLIAN  ELECTRON  DISTRIBUTION  FUNCTIONS  IN 

Z-PINCH  PLASMAS 

1.  Introduction 

A  knowledge  of  the  electron  distribution  function  is  necessary  for  accurate 
determination  of  z-pinch  fluid  and  radiation  dynamics.  Distribution  function  dynamics 
underlie  the  magnetohydrodynamic  behavior  of  the  pinch,  affecting  implosion  time  and 
implosion  energy  as  well  as  the  coupling  efficiency  between  generator  and  load.  The 
radiation  yield  of  the  pinch  may  depend  sensitively  on  these  factors,  which  determine 
the  energy  input  to  the  pinch  and  thus  influence  the  energy  spectrum  of  the  radiation; 
in  addition,  the  flux  of  high-energy  electrons  is  an  immediate  influence  on  the  radiation 
spectrum  due  to  collisional  excitation  in  the  pinch. 

Study  of  the  electron  distribution  function  is  the  only  way  to  investigate  non- 
t hernial  phenomena  like  runaways  and  the  plasma-wave  turbulence  of  anomalous  resistivity. 
Runaway  electrons  have  been  observed  in  z-pinch  experiments,1  despite  the  strong  toroidal 
magnetic  fields  at  the  pinch  edge,  and  may  be  an  important  factor  in  both  fluid  and 
radiation  dynamics,  reducing  coupling  by  providing  an  alternative  current  path  and 
contributing  to  radiative  instabilities,  since  they  are  associated  with  “bright  spots”.1 
But  the  production  and  propagation  of  z-pinch  runaways  is  not  well  understood  at  this 
time.  Anomalous  resistivity  and  the  limiting  of  electron  drift  velocity  due  to  particle- wave 
interactions  play  an  important  role  in  the  pinch,  by  both  increasing  the  rate  of  diffusion  of 
current  and  fields  into  the  plasma  (thereby  altering  the  force  structure  and  so  the  implosion 
dynamics)  and  by  increasing  the  energy  input  to  the  plasma  following  its  assembly  on  axis. 
Information  is  now  needed  on  how  to  control  these  factors  for  appropriate  radiation  yields 
in  a  variety  of  pinch  parameters  and  materials. 

The  shape  of  the  high-energy  distribution  of  electrons  plays  an  important  role  in 
determining  kilovolt  x-ray  conversion  efficiencies  in  z-pinches.  Still,  the  literature  provides 
an  unclear  picture  of  the  effects  that  different  interactions  in  a  z-pinch  have  on  this  part  of 
the  distribution,  since  these  effects  depend  diversely  on  the  state  of  the  plasma  and  on  the 
rate  of  energy  input  or  output  from  it.  Under  various  conditions  the  tail  of  the  distribution 
may  be  enhanced  or  depleted,  may  be  extended  or  may  acquire  a  high  energy  plateau. 
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Most  prior  studies  of  the  Fokker-Planck  equation  have  been  related  to  tokomacs  or 
inertial-confinement  fusion.  In  the  former  case,  the  electric  and  magnetic  fields  are  taken  to 
be  parallel  and  runaway  rates  are  computed  for  those  conditions.  In  the  latter  case,  cross- 
field  transport  is  computed,  but  with  the  assumption  of  weak  magnetic  fields.  Usually,  the 
Fokker-Planck  equation  is  linearized  about  a  Maxwell-Boltzmann  equilibrium  distribution. 
An  earlier  z-pinch  study,  related  to  this  one,  investigated  equilibrium  (time-independent) 
solutions  to  the  nonlinear  Fokker-Planck  equation  in  a  z-pinch.2  These  earlier  approaches 
have  limited  usefulness  for  studying  dynamic  z-pinch  implosions,  which  are  characterized 
by  rapid  change  and  strong,  perpendicular,  spatially  varying  electric  and  magnetic  fields. 

The  principal  effects  that  must  be  accounted  for  in  a  kinetic,  Fokker-Planck  description 
of  the  electron  distribution  function  dynamics  are:  (1)  the  strong  electric  field  heating 
(ohmic  heating);  (2)  the  equally  strong  electron-electron  elastic  collisions;  and  (3) 
the  cliffusion-like  electron-plasma  wave  interactions  that  determine  the  strength  of  the 
instability-driven  plasma  turbulence.  In  this  report  a  basic  theory  is  set  up  which  is 
suitable  for  treating  all  three  of  these  effects  under  the  dynamically  evolving  conditions 
found  in  a  z-pinch  implosion.  However,  only  the  effects  of  ohmic  heating  and  electron- 
electron  collisions  will  be  studied  in  this  report;  the  theory  of  the  turbulence  interaction 
and  its  contribution  to  the  picture  presented  here  will  be  discussed  in  a  later  report.  A 
time-dependent  theory  is  essential  because  the  different  processes,  with  different  timescales, 
change  in  relative  importance  during  the  course  of  the  implosion,  and  because  some  of 
the  phenomena  studied,  like  runaways,  are  essentially  non-equilibrium  in  nature.  Also, 
with  the  time-dependent  theory  it  is  possible  to  continually  assess  the  accuracy  of  the 
approximations  used  on  physical  grounds;  this  provides  a  check  on  the  results. 

The  Maxwell-Boltzmann  equilibrium  distribution  is  produced  by  small-angle  electron- 
electron  collisions,  and  this  equilibrium  will  be  maintained  on  their  (<  0.1  ns)  time 
scale.  Heating  and  cooling  processes  acting  preferentially  in  one  energy  range  will  disrupt 
this  equilibrium.  Because  the  electron  collision  time  decreases  with  velocity,  return  t<> 
equilibrium  is  especially  slow  in  the  distribution  tail.  The  heating  and  cooling  processes 
considered  here  include  ohmic  heating,  compression,  and  radiation.  Ohmic  heating  results 
from  electric  field  acceleration  of  electrons.  Compressional  heating  or  cooling  is  calculated 
from  spatial  gradients  in  the  fluid  velocity,  which  are  supplied  as  input  (e.g.,  from  MHD  z- 
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pinch  simulations).  Radiation  losses  can  be  calculated  from  inelastic  electron-ion  collisions, 
assuming  an  optically  thin  medium. 

This  report  forms  a  compendium  of  analytic  results,  drawn  from  many  sources, 
which  are  relevant  to  z-pinch  kinetics.  It  draws  connections  to  other  work,  particularly 
on  laser  heating  by  inverse  bremsstrahlung  and  on  electron  distributions  in  tokomacs. 
Analytic  results  are  obtained  whenever  possible,  and  individual  effects  on  the  electron 
distribution  are  analyzed  in  depth.  In  particular,  expressions  for  the  conductivity  and 
runaway  production  rates  are  given.  Finally,  a  numerical  solution  is  obtained  to  the  fully 
non-linear  time-dependent  Fokker-Planck  equation.  This  numerical  solution  will  later  be 
used  in  conjunction  with  an  MHD  code  to  estimate  the  magnitude  of  non-thermal  effects 
in  actual  z-pinch  implosions. 


2.  Physical  Picture 


The  system  is  composed  of  electrons,  of  mass  m,  and  various  ion  species  «,  each  with 
mass  Ma  m,  charge  Za  and  density  nQ.  The  mean  charge  number  is  Z  =  ne/n ,  (=s  10), 
where  ne  =  V \QnaZa  and  n;  =  £ lQna.  The  electrons  are  described  by  a  distribution 
function  /(r,u,f).  The  ions  are  assumed  to  maintain  a  Maxwell-Boltzmann  distribution. 
Electrons  are  subject  to  both  electrostatic  forces  mai  =  ZeE  and  magnetic  forces 
muj  =  (Ze/c)B]  the  fields,  E  and  B,  are  locally  uniform  in  time  and  space.  (The  electric 
field  E  should  not  be  confused  with  the  average  electron  energy  {E)). 

The  collision  frequency  for  elastic  Coulomb  scattering  (Rutherford  frequency)  of  a 
particle  of  species  a  moving  with  velocity  vQ  through  a  number  density  rip  of  species  fS 
particles  is 

va0(va)  —  ^n-/3t>a(7rfa  )j  ^(£<*2$)  log  Aa  ,  (1) 


where  Ca  =  2e2/mat>2 ,  called  the  Landau  length,  is  the  distance  of  closest  approach  of  (he 
scattering  particle  and  the  threshold  for  large-angle  collisions;  and  the  Coulomb  logarithm 
is  given  by 


log  Aq  =  log 


3  (fcT)3/2 


(2) 


[2y/ir  Zae3s/nZ 

The  first  bracketed  term  in  eq.  (1)  is  the  cross-sectional  flux  (the  number  of  particles 
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per  second  passing  through  particle  a’s  path),  and  the  second  term  is  the  Coulomb 
scattering  factor,  giving  the  interaction  strength  between  the  two  scattering  particles.* 
Substituting  for  t a  in  eq.  (1)  and  simplifying  gives 

4ffe4  ZlZ%n0  ,  A 

Ua<3~  v3  °gAo'  (3) 

The  “temperature”  of  the  non-equilibrium  system  is  defined  to  be  proportional  to  the 
average  energy  in  the  system: 

kT  =  -(B).  (4) 

9 

Different  proportionality  constants  are  used  by  different  authors;  the  choice  consistent  with 
the  expression  for  temperature  in  an  equilibrium  (3D)  system,  however,  is  q  =  3/2.  The 
thermal  velocity  is  also  defined  in  terms  of  the  average  system  energy: 

-  JM.  (») 

V  m 

This  is  related  to  the  temperature  with  the  constant  q :  if  q  —  3/2,  vth  =  y/okT/ni.  The 
.dectron-electron  collision  frequency  at  the  thermal  velocity  will  be  denoted 

l/R  =  V"(vth).  (6) 

The  geometry  of  the  system  is  that  of  a  z-pinch,  cylindrically  symmetric  about  the 
pinch  axis.  The  system  is  assumed  to  be  spatially  uniform  in  the  <f>  and  £  directions; 
important  radial  gradients  do  exist,  however.  Radial  variations  are  included  in  the  model  by 
dividing  the  pinch  into  concentric,  spatially  uniform  thin  cells,  with  time- varying  inner  and 
outer  radii  and  a  mean  velocity  V(r, t).  The  electron  distribution  function  f{v,t)  is  then 

*  Other  illustrative  forms  for  i/a0  can  be  obtained,  such  as,  for  atomic  physics 
applications, 


=  »,»„<»«; )  (  Jv°J^Ry  )  log  A„, 


where  a0  is  the  Bohr  radius  and  Ry  =  e2 /2ao  is  the  Rydberg  energy;  or,  for  high-energy 
(non-relativistic)  physics,  a0  can  be  replaced  by  rt  =  e2 /me2,  the  classical  electron  radius, 
and  2 Ry  by  mec2,  the  electron  rest-mass  energy. 


to  be  found  within  the  reference  frame  of  each  dynamically  evolving  cell.  This  approach 
includes  the  important  effects  (compressional  heating  or  cooling)  of  bulk  electron  motion 
on  the  electron  dynamics,  without  requiring  that  the  Fokker- Planck  kinetic  equation  be 
solved  in  configuration  space. 

3.  Basic  Equations 

In  this  section  the  model  kinetic  (Boltzmann)  equation  will  be  presented,  an  expansion 
will  be  performed  and  the  first-order  terms  obtained,  and  a  zero-order  nonlinear  Fokker- 
Planck  equation  will  be  obtained  for  analysis. 

Relative  to  a  stationary  observer,  the  electron  distribution  function  satisfies  the 
(collisional)  Boltzmann  equation3 

~  +vL  •  V/  +  (aL  +vL  x  £)■  VVLf  =  .  (7) 

where  vi  is  the  particle  velocity  in  the  stationary  “lab”  frame,  ai  =  ZeE/m  and 
u>  =  (Ze/mc)B.  In  each  z-pinch  cell,  the  system  is  moving  with  a  time-dependent,  average 
fluid  velocity  V[r,t).  The  particle  velocity  and  electric-field  acceleration  in  this  reference 
frame  are 


v(f, 

vi, t)  = 

VL-V(fj) 

(8) 

Ze  -  - 

dV 

a  = 

—  E  +  V  x  ut  - 
m 

'  dt  5 

(9) 

the  comoving  time  derivative 

d/dt 

=  d/dt  +  V  ■  V 

(10) 

gives  changes  with  respect  to  the  moving  frame.  Derivatives  in  the  moving  frame  are  given 

by 

df( v)  df  df  dv  _df  dV^, 

dt  dt  +  dv  dt  df  df  "* 

Vf(r,v,t)  =  V/  —  (W).Vr/.  (11) 

Therefore,  in  a  reference  frame  moving  with  velocity  Vr,  the  Boltzmann  equation  becomes1 
-£  +  v  •  V/  +  [a  +  w  x  w]  •  VJ  -  v  •  W  ■  V„/  =  C'i  +  Cinei  (12) 
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The  last  term  on  the  left-hand  side  vanishes  if  the  mean  velocity  V  is  spatially  uniform; 
it  is  through  this  term  that  compressioual  effects  enter.  Collisions  are  represented  on  the 
right-hand  side  of  the  Boltzmann  equation:  Cti  gives  the  effects  of  elastic  collisions,  a 
Fokker-Planck  term  in  this  model,  and  C\nei  describes  inelastic  processes.  These  terms 
will  be  discussed  in  detail  later. 

Cartesian  tensor  expansion 

Instead  of  being  a  scalar  function  of  the  vector  position  and  velocity  variables,  the 
distribution  function  f(x,v,t)  can  be  expressed  in  terms  of  scalar,  vector,  tensor,  etc. 
functions  of  scalar  position  and  velocity  variables.4  If  only  first-order  terms  are  retained, 
the  expansion  reads 

f{x,v,t)  =5  /o(x,r)  +  ^  v  ■  fi{x, t»)]  (13) 

This  first-order  expansion  is  sufficient  to  determine  averages  of  any  quantity  that  is  a  linear 
function  of  the  scalar  v  and  vector  v.  It  can  always  be  used  if  /i  (and  higher-order  terms) 
is  sufficiently  small,  and  with  certain  geometries  it  can  be  used  regardless  of  the  magnitude 
of  /i .  The  different  terms  of  the  expansion,  when  known,  can  be  used  to  compute  average 
values  such  as  the  density  and  temperature.  For  example,  the  average  value  of  a  function 
<j>{v)  is 

r°° 

(4)  -  —  /  4>fov 2  dv  (14) 

n  Jo 

and  the  average  value  of  a  function  4'(v)  which  is  linearly  dependent  on  the  vector  v  is 

W  =  3n/  4'{f\{v))v2  dv.  (15) 

In  particular,  the  mean  number  density  n(x)  =  (1)  and  current  density  j  =  (env)  are 

n(x)  =  4rr  f  v2f0(x,v)dv.  (16) 

Jo 

Air  f°° 

](*)  =  J  v*fi(*,v)dv.  (17) 

When  the  expansion  of  eq.  (13)  is  substituted  into  the  Boltzmann  equation,  eq.  (12). 
and  the  various  angular  moments  are  taken  (this  procedure  is  illustrated  in  Appendix  II), 
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the  result  is4: 


—y-s-h 

^7  +  rV/0  +  a^+«x/,-  VV  ■  /j 
dt  av 


=  CFP  +  C’ 


(IRa) 


(W  +  V^  +  I3V.V1.^(^)=C1 


186) 


where  Cpp  and  Ci  come  from  Cei,  and  C*  comes  from  C'ine/,  on  the  right-hand-side  of 
eq.  (7).  There  is  no  contribution  to  the  fi  equation  from  inelastic  collisions  because  these 
collisions  are  assumed  to  be  isotropic  and  so  do  no  change  the  average  momentum  of  the 
system. 

The  physical  content  of  eqs.  (18)  is  fairly  clear,  if  one  bears  in  mind  that  fo  describes 
the  particle  density  and  f\  its  current.  Eq.  (18a)  describes  the  change  in  time  of  /<>  resulting 
from  the  net  current  inflow  (the  V  •  /j  and  V  •  V  terms),  work  done  on  the  current,  by  the 
electric  field  (c-/j),  and  collisions,  both  elastic  ( Cpp )  and  inelastic  (C*).  Eq.  (18b)  gives 
changes  in  / 1  corresponding  to  flows  driven  by  density  gradients  (the  V /„  and  VV"  terms) 
and  by  the  electric  field,  which  moves  nonuniformities  in  f(v)  through  velocity  space  (the 
adfo/dv  term).  The  effect  of  these  driving  factors  is  modified  by  the  deflecting  nature  of 
the  magnetic  field  (u;  x  /j),  and  by  elastic  collisions  (Ci). 

Eqs.  (18)  are  rigorous,  except  for  the  neglect  of  higher-order  expansion  terms  like  /; 
in  eq.  (18b).  The  procedure  now  is  to  solve  eq.  (18b)  for  /j  under  physically  reasonable 
simplifying  assumptions,  and  then  to  substitute  the  resulting  expression  for  a  •  /j  into 
eq.  (18a)  to  get  an  equation  that  can  be  numerically  solved  for  /0.  An  immediate 
simplification  comes  if  we  ignore  the  (assumed)  weak  effects  of  the  currents  that  are  driven 
by  density  and  temperature  gradients,  i.e.,  only  gradients  in  the  fluid  velocity  V  will  be 
retained.  Then  V  •  fi  =0  and  V/0  =  0,  and  eqs.  (18)  may  be  summarized  as  follows: 


dfo 

dt 


=  CFP  +  C"  +  C'c  +  Ca 


l  —  ^  f 0  ,  -*  p  rp  ,  /T 

+  a  a  — k  w  x  /j  —  T  +  C  i 
dt  a  v 


(19r») 


(106) 


where  Cpp  describes  elastic  collisions,  C*  describes  inelastic  collisions,  C’c  and  T  contain 
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compressional  effects,  and  Ca  is  due  to  electric-field  heating.  These  terms  are  discussed  in 
detail  below. 


Elastic  collisions 

The  elastic  collision  terms  in  eqs.  (19)  come  from  the  expansion  of  the  Cei  term  in 
eq.  (12),  which  includes  contributions  from  all  particle  species:  C'i  =  'LQC\V-  The  zero- 
order  term  Cpp  describes  energy  exchange  in  collisions,  while  the  first-order  term  describes 
momentum  transfer,  and  so  has  no  contributions  from  like-particle  collisions.  If  FQ,  is  the 
i-th  order  term  in  the  expansion  of  the  distribution  function  for  species  a,  the  general 
expressions  are4,5  [The  Cpp  equation  is  derived  in  Appendix  I]: 

p(“)  _  &  rn  o  ,  v  dfo  (  t o  ,  to  \1  ,9n  » 

“  V2  to  [a4/o/°  +  2  (20fl) 

^(o)  _  lea  f  t'~  9‘fi  0  0  V  8f±  /3rn  „  rO,9Ao\,1//  orO  ,  f0  9r0  , 

Cl  “^\3  dv'{ h  +  -l)+  3^7  +2J-')  +  3/l(_3/° 


+  47r1’3  — -(/iFqo  +  Faifo)  -I-  —  -4°  (Q  +  J- 2) 
M  5  ov~ 


(206) 


with  the  definitions 


Yq0  =  4tt  ( ZQZ0e2/mo )2  log  Aa  = 

n/J 


i]i(v)=^L  F°'(v)vj+2dv 

Jjl(v)= ^7  fv°°  f*'{v)vj+2  dv 


These  expressions  can  be  simplified. 

Electron-electron  collisions  dominate  in  Cfp  because  of  the  extreme  mass  difference 
between  electrons  and  ions  (just  as  there  is  very  little  collisional  interaction  between  a 
basketball  and  a  Mack  truck,  which  have  about  the  same  mass  ratio).  The  electron-ion 
term  Cpp)  is  proportional  to  m/il/,  and  so  is  ignored.  With  only  electron-electron  collisions 
included,  the  term  Cpp  in  eq.  (20a)  becomes 

CFP  =  47r^T^;  |/o  J  fo(v')v'2dv '  +  (~2  fo  dv' +  v  J  f0(  (’>'  rfr'j 
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In  the  Ci  term  it  is  ion-electron  collisions  that  are  important,  because  no  momentum 
is  exchanged  in  collisions  between  like  particles.  In  eq.  (20b),  the  large  ion  mass  allows 
us  to  approximate  Fn  %  0  -  the  ions  acquire  no  motion  from  electron-ion  collisions.  If 
the  ions  are  assumed  to  be  in  a  spherically  symmetric  Maxwellian  distribution,  then  C\ 
simplifies  considerably.  The  integrals  involved  in  I ■  and  Jj  then  become  functions  of  the 
error  function  $(r y/ M /2kT)  and  its  derivatives4,  and  because  of  the  relatively  large  ion 
mass  only  the  large-argument  limit  of  these  functions  need  be  considered.  In  this  limit, 
Jl\  ~  0;  /0°  ~  rii\  and  i2°  %  3-riikT /Mv2 ,  and  if  terms  proportional  to  m/M  are  again 
dropped,  eq  (20b)  becomes  (see  eq.  (3)) 

C,  =  -vjx  (25) 

Essentially,  these  approximations  assume  that  the  ions  affect  but  are  not  affected  by 
the  electrons.  Ion  equilibrium  is  plausible  because  of  the  large  ion-ion  collision  rate 
(un/ute  s;  Z4)  as  well  as  because  of  the  large  ion  mass. 

Because  electron-ion  collisions  are  much  more  frequent  than  electron-electron  collisions 
(V'i/vee  =  Z ),  and  because  f\  changes  on  the  timescale  of  vcl  while  /0  changes  like  vee, 
then  fi  equilibrates  much  faster  than  /o.  This  ordering  of  time  development  will  be  used 
later.  The  electron-ion  collisions  represented  by  uei  damp  any  increase  in  fi  and  cause  a 
drag  on  electron  flow;  the  decrease  of  with  increasing  v  results  in  high-energy  electrons 
becoming  “runaways”. 

Inelastic  collisions 

The  inelastic  collision  term  C*  in  eq.  (19a)  can  be  obtained  by  requiring  conservation 
of  particle  number  in  collisions.  Inelastic  collisions  transfer  energy  from  electrons  to  ions, 
exciting  the  ions  to  a  higher  energy  level.  The  ions  are  assumed  to  immediately  re-radiate 
this  energy  through  the  optically  thir  medium,  so  that  electron-ion  collisions  alone  do  not 
significantly  change  the  ion  level  populations.  The  minimum  free-electron  energy  to  excite 
a  transition  from  level  a  to  level  b  is  the  excitation  energy  The  number  of  ion  levels 
contributing  significantly  to  the  electron  kinetics  is  generally  small,  of  the  order  of  10. 

Since  ionization  and  recombination  are  not  considered  here,  the  number  of  ions  and 
free  electrons  is  conserved  in  inelastic  collisions.  C'ollisional  losses  are  described,  then,  by 


9 


subtracting  electrons  from  a  higher  energy  region  of  the  distribution  function  and  replacing 
them  at  a  lower  energy  region.  Each  of  the  possible  inelastic  scattering  processes  can  be 
numbered  by  an  index  j\  if  the  collision  frequency  for  process  j  is  v*-  ,  then  the  collision 
term  is  of  the  form 

C%  =  ^  [_^*(f)/o(€)  +  \j — t/J(e  +  (j)fo(e  +  Cj)]  (26) 

j 


The  first  term  in  the  sum  represents  the  rate  of  scattering  to  lower  energy  of  particles 
with  energy  e,  while  the  second  term  represents  the  influx  of  particles  to  energy  c  from 
collisions  involving  higher-energy  particles.  C*(<)  thus  gives  the  net  rate  of  increase  in 
the  distribution  function  at  energy  e  due  to  all  inelastic  scattering  processes.  The  factor 
y/(t.  +  tj)/ 1  makes  allowance  for  the  greater  phase-space  volume  available  at  the  energy 
6  +  Cj  from  which  the  electrons  scatter  to  energy  c:  the  number  of  electrons  scattered 
between  energy  e  and  e  4-  ej  is  always  v* (e)fo{t)  4ir y/e  de. 

The  inelastic  collision  frequency  for  an  electron  (velocity  v)  to  excite  an  ion  from  level 
a  (density  /ant)  to  level  b  is 


Kb  =  /«n« 


flqb(  t)/ 9a 

KRy 


* 

0(e  -  tab) 


t>. 


(27) 


The  quantity  in  brackets  above  is  the  inelastic  cross-section,  which  is  expressed  in  terms  of 
the  measured  dimensionless  collision  strength  fl0j,(e).  The  multiplicity  of  level  a  is  ga.  The 
Heaviside  step  function  6  relects  the  energy  threshold  for  the  transition.  A  plot  of  typical 
collision  strengths  for  argon  is  given  in  fig.  (1);  typically,  Uai,  is  of  order  1,  varying  little 
over  the  energy  range.  The  ratio  of  the  inelastic  collision  frequency  to  the  electron-electron 
collision  frequency  uee  is 


uab  __ 

V" 


fa  flat 
9a 


( 


*/Ry 

4  Z  log  A 


ea6  )• 


(28) 


This  is  generally  small;  since  the  only  terms  not  of  order  one  are  Z  and  log  A,  which  are 
each  near  10,  an  effect  of  perhaps  a  few  percent  for  each  transition  might  be  anticipated 
from  inelastic  collisions.  Because  the  ratio  increases  with  e,  however,  greater  effects  should 
appear  in  the  tail  of  the  distribution,  in  particular  at  energies  greater  than  Z  times  the 
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thermal  energy.  Also,  since  the  inelastic  term  is  a  sum,  the  cumulative  effect  of  inelastic 
collisions  could  be  quite  significant.,  if  many  important  transitions  are  available  to  the 
system. 


Compression aJ  heating 

The  compressional  term  in  eq.  (19a)  comes  directly  from  eq.  (18a): 

-(*•*)§&  (29» 

Tliis  describes  changes  in  the  distribution  function  fo  caused  by  a  plasma  compression  or 
expansion  which  is  reflected  in  the  nonuniform  fluid  velocity  V .  Taking  moments  of  Cc 
reveals  that  this  term  conserves  neither  local  particle  density  nor  average  energy: 

,301 

(*!r)c --<*>!*•*  1311 

The  change  in  density  in  eq.  (30)  is  just  that  resulting  from  the  continuity  equation  for 
compression  or  expansion  of  a  uniform-density  flow  of  electrons.  The  change  in  energy 
density  in  eq.  (31)  is  caused  by  PdV  work:  the  pressure  P  =  2(E) /3,  and  integrating  V  •  V' 
shows  it  to  be  the  time  rate  of  change  of  the  system  volume. 

The  effect  of  compression  on  the  first-order  equation,  eq.  (19b),  conies  from 

f  =  VV-h  +[W  +  Wt  +  I2VT].  ^--1-  (  ]  ,  (32) 

5  ov  \  v  J 

where  I2  denotes  the  3x3  unit  tensor.  This  term  is  not  simple;  in  a  cylindrical  coordinate 
system  and  with  V  in  the  r  direction  only,  the  components  of  T  are 


Since  there  is  no  shear  in  the  cylindrically-symmetric  flow,  no  mixing  of  the  different 
components  of  fi  occurs  as  a  result  of  the  compression.  The  effects  included  in  T  are 
small,  since  they  are  second  order  in  the  small  quantities  f\  and  VV.  They  modify  the 
response  of  the  fluid  to  the  applied  electric  field  due  to  a  changing  plasma  density,  and 
should  be  small  compared  to  t/ej.  In  this  treatment,  therefore,  we  take  T  =  0. 


Electric  field  -  Ohmic  Heating 

Work  done  on  the  electrons  by  the  electric  field  shows  up  in  the  zero-order  distribution 
function  through  the  term  Ca  in  eq.  (19a): 


Ca 


1  d 


3i>2  dv 


-  v-a-  /l]  . 


(33) 


When  this  energy  is  distributed  throughout  the  plasma  by  electron-ion  collisions,  it  causes 
the  bulk  heating  which  is  ohmic  heating.  The  magnetic  field  is  important  in  ohmic  heating, 
since  is  alters  the  electron  trajectory.  It  will  be  seen  that  this  shifts  the  heating  to  lower- 
energy  electrons;  the  heat  then  must  be  transferred  to  high-energy  electrons  by  elastic 
scattering,  whose  strength  decreases  with  energy. 

To  find  C a ,  eq.  (19b)  must  be  solved  for  a  ■  f\.  This  is  straightforward  in  the  center- 
of-mass  system.  With  eq.  (25)  for  C\  and  T  =  0,  eq.  (19b)  is 

-£■  +<Jj  x  £  =  -  a~  (34) 

An  equilibrium  solution  is  sought  for  f\  {fi  reaches  equilibrium  on  the  timescale  of 
i/ci,  which  is  Z  times  faster  than  the  ute  timescale  of  /0).  It  can  readily  be  verified  that  if 
dfi  / dt  =  0,  the  equilibrium  /j  is 


/i  = 


w2  +  v-- 


_  _  u>  •  a  _ 

u >  x  a - 1 v  —  u.,a 


dfo 

dv 


Thus 


a  •  /i  =  - 


a2  +  (a  •  v /vei)2  df0 


(35) 


(36) 


i^(l  +ur/i/;.)  dv 

When  the  electric  field  acts  parallel  to  the  magnetic  field,  C'a  is  independent  of  u>,  as 
it  should  be.  When  the  electric  and  magnetic  forces  are  perpendicular,  as  in  the  Bennett 
z-pinch,  then  a  •  u?  =  0,  and  the  olunic  heating  term  Ca  becomes 

a2  dfo  1 


^  i  a 

c„  =  cEB  =  — - 


i/e,  (1  4-  u>2/i/2t)  dv  J 


(37) 
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Normalization 

For  ease  of  computation,  it  is  useful  to  express  the  quantities  given  above  in  terms 
of  dimensionless  variables,  using  characteristic  velocity,  time  and  space  variables  for 
normalization.  A  dimensionless  velocity  representation  is  presented  and  then  used  as  a 
stepping  stone  to  the  dimensionless  energy  representation  used  in  the  remainder  of  this 
work. 

Characteristic  magnitudes  used  in  scaling  are:  for  time,  the  thermal  electron-electron 
collision  time  tr  =  i/^1  (see  eq.  (6));  for  length,  the  mean  electron  separation  n-1/3;  for 
velocity  or  energy,  the  thermal  velocity  v1h  (see  eq.  (5)).  The  characteristic  electric  field 
is  the  Dreicer  field10  Ed  =  which  is  the  electric  field  that  would  accelerate  an 

electron  from  rest  to  the  thermal  velocity  in  one  thermal  collision  time.  The  characteristic 
magnetic  field  has  a  cyclotron  frequency  of  i/r:  Bd  —  mcu^/e. 

The  dimensionless  velocity  £  is  defined  with  respect  .to  the  thermal  velocity  (eq.  (5)): 

(  =  —■  (38) 

vth 

It  can  be  seen  from  eq.  (16)  that  f0(v)  has  dimensions  of  •  Therefore  the 

dimensionless  distribution  function  /0  is: 

v3 

Jo  =  —fo{  v)  (39) 

71 

The  function  /o  is  a  probability  distribution  whose  average  density  (eq.  (16))  is  normalized 
to  one.  To  make  eq.  (19a)  dimensionless,  both  sides  are  multiplied  by  v^tr/t?.  The  new 
elastic  collision  term  is  then 

Cfp  =  X~i~-CFp  =  ’  (4°) 


where 


The  dimensionless  inelastic-collision  term  C *  is  just  eq.  (26)  with  the  collision 
frequency  Uj  replaced  by  the  dimensionless  frequency 


(43) 


Since  the  compression  term  Cc  does  not  conserve  density,  a  term  proportional  to 
dn/dt  results  on  substituting  the  dimensionless  distribution  function  /o  into  eq.  (19a).  If 
this  term  is  evaluated  from  eq.  (30),  it  can  be  combined  with  Cc  and  the  normalized 
kinetic  equation  will  then  maintain  a  constant  normalization.  The  true  density  must  then 
be  computed  from  the  integral  of  eq.  (30).  The  dimensionless  compression  term  is  then 


Cc  ~ 


vthTR  c  1  7  _  ^  V 13-  ^  /v3  (  \ 


n 


H2  d( 


The  dimensionless  ohmic  heating  term  is 


C  eb  = 


dCU) 

C-  d( 


where 


(44) 


(45) 


q  =  ze- 

i 


C(0  = 


3Z2  i  +  w^/z2' 


(46) 

(47) 


and  the  dimensionless  electric  and  magnetic  field  variables  are  £  =  Ef  Ed  and  Q  =  ui / ur. 
The  velocity  dependence  of  the  collision  frequency  uei(v)  has  been  used: 


—  =  n1n.r,  =  zf-3 

i/fl  nelee 

In  terms  of  dimensionless  velocity  variables,  the  kinetic  equation,  eq.  (19a)  is 


(48) 


Hi 

dr 


=  1— ( 


N(t)+  yrRV-r 


/o(0 +  [/>(€) +  aC(0] 


+  c\ 


19) 


The  dimensionless  energy  variable  is  defined  in  terms  of  the  temperature  (see  eq.  (4)): 

m  i'2/2 


€  = 


kT 


=  qE 


50) 
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Important  relations  between  the  dimensionless  velocity  and  energy  variables  are 

^  =  d(  =  2^fC-  (!il) 

In  terms  of  the  dimensionless  energy,  and  with  the  following  definitions  (new  definitions 
from  €-dependent  functions  are  denoted  by  a  new  argument): 

N(e)  =  4tt  f  de' Vd  f0(e' )  (52) 

Jo 

V(e)  =  4tr |  \j  +  <3/!  J  dt'fo(t')  153) 

a  =  Z£’  =  -  -ff  ,  E’  (54o) 

2m(E)Ufl 

((e>  =  az-^i  +  nv/zv  (r>'1(,) 

C-(€)=  +  |  (55) 

a  states  l  \  e  l /R  J 

"aVO  g3/2«rs  f/.n«fc(0l  ir  L, 

^a6(«)  =  -  =  7Ty\ - 7  -  - T~r  °(e  ~  fab)  (56) 

ur  2ZlogA  g  a  mc-y/e 


(546) 


;C  +  ^a%(e  +  fab) 


/o(e  4-  €„&) 


.  <■  /  v  1  abV1  /  1  FS  ja*“ao\'-r  \  ai  \  <rn 

^ab(e)  =  -  =  777 - 7  -  - —r =  0«-e.b  (56) 

UR  2Z  log  A  [  pa  [mc-y/cj 

where  aps  ~  1/137  is  the  fine-structure  constant,  the  kinetic  equation,  eq.  (19a)  becomes: 


dfn  1  d  ff 

77  =  -^  ^(0  + 


2tr  d(rV 
3r  Or 


-<J,!  /.(0  +  [i>(«)  +  oi(0!^}  +  c' 


dr  So  l  "  3r  dr‘  j  J 

This  form  of  the  Fokker-Planck  equation  will  be  the  basis  for  the  following  analysis 


(<)  (57) 


4.  Analytical  Results 

In  this  section,  analytical  expressions  are  found  for  the  conductivity  both  parallel  to 
and  perpendicular  to  the  magnetic  field;  different  expressions  are  presented  for  studying 
“runaway”  in  a  magnetic  field,  self-similar  solutions  are  found  for  the  time-development 
of  the  distribution  function,  and  the  effects  of  ohmic  heating,  compression,  and  inelastic 
collisions  are  studied. 

Conductivity 

The  conductivity  of  the  plasma  can  be  computed  by  using  the  definition  of  the 
conductivity  tensor  er: 

J  =  cr  ■  E  (58) 
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with  the  current  j  computed  from  the  distribution  function  according  to  eq.  (17). 

With  cylindrical  symmetry  and  a  nonzero  magnetic  field,  the  conductivity  tensor 
has  two  different  components,  a  parallel  conductivity  <t\\  along  the  electric  field,  and  a 
perpendicular  (or  Hall)  conductivity  <rj_  perpendicular  to  the  electric  field.  The  parallel 
conductivity  is  easily  obtained  from  eq.  (36): 


ill  = 


E  r  <■*  sf» . 

3 mniYei  J0  1  +  u6u ;2/n?Ke2  dv 


(59)) 


Identifying  terms,  and  expressing  the  integral  in  the  dimensionless  energy  form,  the  parallel 
conductivity  is 

ZvR(»,\*r  e3  dfp 

\u >)  J0  Z2q3/U2  +  e3  de  d  (60 


<T||  -  -■ 


where  u>p  =  y/Arrne2 /m  is  the  electron  plasma  frequency,  and  fl  =  u/vr,  as  defined  in  the 
text.  The  Hall  conductivity  can  be  obtained  directly  from  eq.  (34),  from  which  follows 
(with  a  ■  uj  =  0) 

ax  fi  =  (“  •  /i)  (61) 

This  gives  the  perpendicular  distribution  in  terms  of  the  parallel  one,  so  that 


<7-L  = 


1  '°° 


/ 


f9/2 


dfo 

3 q3/2  «  J0  Z2q3/Q2  +  c3  de 


de. 


(62) 


Both  the  magnetic  field  and  the  temperature  affect  the  conductivity,  as  does  the  electric 
field  because  of  its  effect  on  the  distribution  function.  In  the  strong  magnetic  field  limit 
(ft  >  Zq3/2)  the  conductivity  reduces  to 


<T||(u  —  oo)  a  2^2  (^)  /o(e  =  0)  — .  0  (63a) 

-  CO)  « -0.  (636) 

In  the  weak  magnetic  field  limit,  the  conductivity  becomes 

w2  r°° 

(7||(u;  —  0)  as  -  p  /  e2/0(e)de  (64n) 

T Z VR  J0 

g~x(<^  — >  0)  a:  u  ^^2g9/2  (~^)  Jo  €‘ /2  fo(()  d*  — *■  0.  (64b) 
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These  agree,  in  the  appropriate  limit,  with  accepted  expressions  for  the  static  magneto¬ 
conductivity9. 

At  low  temperatures,  the  runaway  effect  dominates,  and  raising  the  temperature 
decreases  collisions  and  increases  conductivity.  At  high  temperatures,  however,  the 
magnetic  field  dominates,  and  collisions  increase  conductivity  by  breaking  the  magnetic 
confinement;  thus,  increasing  the  temperature  and  decreasing  collisions  actually  decreases 
the  conductivity.  This  behavior  is  shown  in  fig.  (2),  where  <T\\  (calculated  with  /0(c)  a 
Maxwellian)  is  plotted  as  a  function  of  temperature. 

The  current  in  the  pinch  is  not  parallel  to  the  electric  field.  Instead,  from  eq.  (58), 

j  =  [f<x j.  +  £cr||]  £,  (65) 

which  is  a  conical  flow  of  current  at  an  angle  9j  =  tan ~1(<Tn/<r±)  from  the  pinch  axis.  In 
the  actual  pinch,  this  non-axial  current  will  probably  be  of  short  duration;  as  with  the 
Hall  effect  in  a  conductor,  a  radial  charge-separation  electric  field  will  quickly  develop  in 
the  pinch  which  will  direct  the  current  back  onto  the  axis.  The  true  parallel  conductivity 
is  here  taken  to  be  (T\\ . 

As  an  alternative  to  the  above  calculation,  the  parallel  conductivity  in  eq.  (60)  could 
also  have  been  obtained  by  computing  the  heating  rate  due  to  Ceb •  The  ohmic  heating  is 
just  equal  to  <T||E2,  and  the  coefficient  of  E2  in  the  resulting  expression  is  just  <T||.  These 
two  methods  are  equivalent  and  self-consistent. 

Runaway  electrons 

The  Coulomb  collision  rate  decreases  with  relative  energy  (see  eq.  (3));  this  implies 
that  particles  with  large  enough  kinetic  energy  are  not  significantly  affected  by  collisions. 
Electrons  with  sufficient  energy  become  runaways,  and  behave  almost  as  free  particles 
accelerated  by  the  electric  field.  With  no  magnetic  field,  runaway  electrons  are  those  with 
a  vel  city  greater  than  the  critical  velocity10  rc  r ihy/oED/qE.  (Dreicer’s  “critical  field" 
Ec  is  Ep/q  in  the  present  notation).  The  critical  velocity  vc  is  independent  of  vih  since 
ur  <x  v^3);  as  the  system  is  heated  and  vth  increases,  an  increasing  number  of  electrons 
have  velocities  above  vc  and  thus  run  away. 
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When  inelastic  collisions  and  a  magnetic  field  are  present,  the  situation  is  more 
complicated.  Inelastic  collisions  may  inhibit  runaways  by  providing  an  energy  sink 
and  so  preventing  sufficient  acceleration.  A  strong  perpendicular  magnetic  field  acts  to 
stop  transport  of  charged  particles  on  any  length  scale  larger  than  the  cyclotron  radius 
rc  =  2.3845 \/K / B  cm  (where  K  is  the  electron  energy  in  eV  and  B  is  the  magnetic  field 
in  Gauss).  If  this  distance  is  smaller  than  the  anode-cathode  gap,  the  electrons  can  only 
flow  from  anode  to  cathode  if  collisions  break  the  magnetic  field  localization.  Thus,  for 
strong  magnetic  fields,  energetic  electrons  are  more  efficiently  trapped  by  the  field  and 
freely-accelerating  runaways  are  eliminated.  Still,  even  in  this  case  energetic  electrons  can 
play  an  important  role  in  the  system,  and  can  acquire  significant  directionality  from  drifts 
or  collisions. 

It  is  clear  that  with  the  z-pinch,  unlike  with  the  tokoinac,  one  cannot  define  an 
arbitrary  energy  cutoff  beyond  which  electrons  behave  as  if  freely  accelerated  by  the  electric 
field.  The  appropriate  definition  of  “runaway”  depends  on  the  application.  Finite-size 
effects  can  be  estimated  using  the  mean  free  path  and  the  physical  dimensions  of  the 
system.  The  rate  at  which  high-energy  electrons  are  produced  can  be  found  by  defining  a 
velocity  threshold  (e.g.,  vc)  and  calculating  the  electron  flux  past  this  threshold.  Beams 
of  electrons  can  be  found  by  looking  for  peaks  in  the  distribution  function.  Because  of  the 
bending  effect  of  the  magnetic  field,  peaks  in  the  distribution  function  rarely  correspond 
to  monodirectional  beams;  there  is  usually  comparable  electron  flux  both  parallel  and 
antiparallel  to  the  electric  field.  Directionality  can  be  checked  using  the  first  term  in  the 
expansion:  /o(t’)  4-  fi(v)  is  the  directional  distribution  function. 

Anisotropy  in  the  electron  velocity  distribution  shows  up  in  fi ,  which  can  be 
considered  a  “current  distribution  function"  (see  eq.  (17)).  A  nonzero  /i(e)  indicates 
a  directed  streaming  at  energy  e.  In  the  presence  of  a  magnetic  field,  /i  is  not  parallel  to 
the  electric  field  a,  though  it  is  always  perpendicular  to  the  magnetic  field  2  (this  follows 
from  eq.  (34),  since  a  -2  =  0).  The  magnetic  field  induces  an  E  x  B  drift  in  the  electron 
motion,  while  electron-ion  collisions  interrupt  the  drift  motion  and  allow  the  current  to 
flow  along  E.  In  a  z-pinch  geometry  2  in  the  <i>  direction  and  E  and  /j  in  the  f-i  plane. 
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the  angle  9(e)  between  /i(e)  and  a  is  given  by 


(66) 


In  the  real  pinch  the  electric  field  is  composed  of  an  applied  field  in  the  z  direction,  and  an 


induced,  r-dependent  component  in  the  radial  direction  to  inhibit  charge  collection  on  the 
edge  of  the  pinch  (see  the  discussion  following  eq.  (65)).  The  flux  of  electrons  at  a  given 
energy,  fi(e),  will  also  generally  have  a  radial  component,  which  tends  to  prevent  focused 
beams  and,  if  the  magnetic  field  is  strong  enough  that  collisions!  electrons  flow  radially, 
to  promote  radial  isotropy  in  the  plasma.  Electron-ion  collisions,  by  interfering  with  this 
magnetic  field-induced  mixing,  tend  to  promote  beam-type  anisotropy  -  just  the  opposite 
of  the  situation  with  no  magnetic  field. 

The  magnitude  of  the  electron  directionality  can  be  found  from  the  full  distribution 
function.  Parallel  to  the  electric  field,  the  distribution  function  is  /0(e)  -f  a  ■  fi(e),  while 
antiparallel  to  the  field  it  is  /0(c)  —  a  •  fi(e).  The  dimensional  expression  for  a  ■  fi(e)  is 


(see  eq.  (36)): 


/!(<)  =  -  ~e 


1  +  Q2e3/(Z2q3) 


Thus,  for  small  magnetic  fields,  f\  oc  e2df0/de,  and  the  distribution  function  becomes 
more  directional  at  higher  energy  (the  runaway  effect),  but  for  large  magnetic  fields, 
f\  oc  fi-2 e-1  dfo /de,  and  the  distribution  function  becomes  less  directional  at  higher 
energy.  The  anisotropic  collision-dominated  low-energy  behavior  changes  to  the  isotropic 
magnetic-dominated  high-energy  behavior  at  the  maximum  of  a  •  /i(<)//MB(f)i  which  for 
a  Maxwellian  /o(c)  is  given  by 

2Z2  1//3 

(max  ~  9  (68) 


Practically,  of  course,  the  magnetic  field  confines  particles  only  within  their  cyclotron 
radius  rc;  whether  any  real  magnetic  field-induced  isotropy  is  obtained,  or  if  instead  a 
significant  loss  of  particles  occurs,  is  determined  by  the  dimensions  of  the  system.  In  a 


19 


finite-size  system,  the  mean  free  path 

p(vth)  = 


Vtk 


(69) 


uR  +  u* 

is  an  important  quantity.  (In  this  equation,  P*  denotes  the  mean  inelastic  collision 
frequency).  When  an  electron’s  mean  free  path  exceeds  the  dimensions  of  the  system, 
the  electron  can  be  said  to  have  “run  away”,  and  the  distribution  function  can  be  cut  off 
at  that  energy.  To  estimate  the  number  of  electrons  whose  free  path  exceeds  the  given 
distance  p ,  eq.  (16)  can  be  generalized,  to  give  a  “runaway”  electron  density  of 


r° o 

np-  4tt  /  dt  \ft  f0(e) 

Jt(p) 


(70) 


where  t(p)  comes  from  inverting  eq.  (69);  if  only  electron-electron  collisions  are  considered, 
then 

«(p)  =  (”i) 

Once  an  energy  is  somehow  chosen  at  which  electrons  are  deemed  to  be  lost  from 
the  system  by  runaway,  the  rate  at  which  electrons  become  runaways  can  be  estimated 
by  calculating  the  flux  of  electrons  through  a  suitable  surface  in  velocity  space.11,12  This 
flux  is  caused  both  by  collisional  diffusion  and  by  the  overall  heating  of  the  system.  The 
runaway  rate  through  a  sphere  in  velocity  space  with  radius  corresponding  to  energy  t  is 

(72) 


rW  = 


n(e)  dt 


where 


n(t)  =  N(t)  =  An  j  dt  V7'  fo(t').  (73) 

Jo 

Note  that  the  total  density  n  =  n(oo)  =  1  in  the  normalized  units  adopted  here.  The 
runaway  rate  can  be  found  by  integrating  from  eq.  (57),  with  N( 0)  =  D(0)  =  C(0)  =  0: 


=  - 


47T 

W) 


/o(e)  + 


+  f  dt  sft'  v*  {t  )fo(t  )\ 


D(t)  +  q((c) 


dfo(t) 

de 


(74) 


For  a  Maxwellian  distribution  the  Fokker-Planck  terms  N(t)fo(t)  and  £>(e)/o  (f)  m  this 
expression  cancel,  and  electron  diffusion  in  velocity  space  is  due  only  to  compressional. 
ohmic  and  inelastic  collision  processes. 
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Self- Similar  Solutions 

The  “standard”  Fokker-Planck  equation  (which  includes  C'pp  only)  causes  the 
electrons  to  approach  a  Maxwell- Boltzmann  distribution  from  any  initial  state.  Those 
Fokker-Planck  terms  are  present  in  the  description  given  here;  acting  against  the  tendency 
to  approach  a  Maxwellian,  however,  are  the  two  new  terms:  the  heating  term  represented 
by  o(,‘(e),  and  the  radiative  cooling  term  represented  by  C*(e). 

The  ohmic  heating  term  (eqs.  (54))  can  be  written  as 

C«(e)  =  AZq3'-g{t)  (75) 


where  A  is  proportional  tow  2 ,  and  the  shape  function 

,3 


9(c)  = 


zv/n2  +  <3 


(76) 


determines  where  the  heating  is  focused.  A  plot  of  <;(e)  and  its  derivative  is  given  in 
fig.  (3).  The  function  starts  at  zero,  rises  abruptly,  and  then  asymptotically  approaches 
its  maximum  value  of  1  at  high  energies.  Low-energy  particles  experience  little  heating, 
because  the  large  low-energy  collision  rates  effectively  randomize  the  electron  motion  and 
prevent  the  electric  field  from  doing  net  work.  High-energy  particles  are  easier  to  heat 
up  because  more  phase  space  is  available  to  them  (this  is  reflected  in  the  factor  of  v2  in 
eq.  (37)). 

Ohmic  heating  is  focused  at  the  inflection  point  of  g(e),  located  at  ema*  (see  eq.  (68)), 
at  the  transition  point  between  the  collisional  and  magnetized  regions  of  the  plasma.  For 
“strong”  magnetic  fields  (f l/Z  >  2),  g{l)  %  1,  the  function  g(e)  is  roughly  constant  over 
most  of  the  distribution,  and  most  of  the  ohmic  heating  occurs  below  the  thermal  energy. 
For  weak  fields,  on  the  other  hand,  the  heating  varies  greatly  over  the  distribution  and 
is  strongest-in  the  tail.  Both  strong  and  weak  field  heating  can  occur  in  a  z-pinch:  the 
approximate  magnitude  of  ft  is 


When  high-magnetic-field  ohmic  heating  dominates  and  p(e)  is  effectively  constant, 
the  distribution  function  has  a  self-similar  form  analogous  to  that  which  has  been  found 
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for  laser-heated  plasmas.6  [See  Appendix  III  for  a  derivation  of  the  laser-heating  inverse 
bremsstrahlung  term.)  When  inelastic  collisions  are  ignored,  the  self-similar  form  is6,7 


/ss(M)  = 


(78) 


y(t)3/2 

where  m  is  determined  by  the  relative  strengths  of  the  heating  and  elastic  collision  terms, 
and  the  normalization  and  kinetic  energy  integrals  determine 

m 


a  = 


y  =  qkT 


8flT(3/m) 
r(3/m) 


(79) 

(80) 


r(5/m) 

The  function  fss  is  an  exact  solution  to  eq.  (57)  with  no  inelastic  or  compressional  terms 
(C*  —  V-Vr  =  0),  in  two  special  cases:  When  elastic  collision  dominate  (a  =  0),  the  solution 
is  a  Maxwell- Boltzmann  distribution,  which  is  just  fss  with  m  =  2  and  y  independent  of 
time;  and  when  ohmic  heating  dominates  (a  —  oo),  the  solution  is  fss  with  m  =  5.  In 
other  cases,  fss  is  an  approximate  solution  to  eq.  (57). 

When  only  ohmic  heating  influences  the  distribution  (a  — ►  oo),  it  is  possible  to 
analytically  determine  the  value  of  m  and  y(t)  by  substituting  eq.  (78)  into  eq.  (57). 
The  result  is: 


dy 

dr 


lysft 


m/2-2 

'm-2  —  m(e/y)m/2' 

\y) 

3  —  m(e/y)m/2 

where 


4q5/2 

A’  =“«»)= 


This  has  a  self-consistent  (i.e.,  e-independent)  solution  only  for  m  =  5: 

2/S 
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y(t)  =  [—Atj  +  const. 


(81) 


(82) 


(83) 


An  analytic  solution  for  m  and  y(t)  is  no  longer  possible  when  a  is  finite,  since  the  elastic 
collision  terms  make  the  equation  nonlinear  and  the  inelastic  collision  terms  are  non-local. 
However,  in  the  limit  a  —  0,  the  m  =  2  Maxwell- Boltzmann  solution  must  be  recovered. 
Numerical  solutions  have  given  a  reasonably  accurate  formula  for  the  exponent  m  in  (lie 
absence  of  inelastic  scattering:' 

3 


m(a)  =  2  + 


1  +  1.66/a0-724 


(84] 
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This  is  plotted  as  a  function  of  the  electric  field  (in  V/cm)  in  fig.  (4). 

The  time  dependence  of  the  kinetic  energy  of  a  self-similar  solution  can  be  found  from 
the  ohmic  heating,  compressional  and  inelasti'  collision  terms  alone,  since  the  Fokker- 
Planck  terms  conserve  energy: 


d(E) 

dr 


=  — 4 tto  f°°  C(e)^r  dt  —  — V  •  V (E)  4-  4rr  re3/-C'(e)de 
Jo  3  Jo 

'  '  a  »«atea  '/0 


=  Air  A  i 


( r~r 

3 


vab(()f(e)y/e  de.  (85) 


where  the  step  function  in  C*  allows  the  lower  limit  of  the  integral  to  be  zero. 

It  is  an  interesting  fact  that  the  self-similar  distribution  is  left  unchanged  by 
compressional  heating  or  cooling.  This  may  be  readily  verified  by  substituting  eq.  (78) 
into  the  density-conserving  compressional  term  from  eq.  (57) 

(el,2/('))  •  ,86) 

where  A2  =  [2rfl/3r]<?(rVr  )/dr.  The  time-development  equation  df/dr  =  Cc  is  then 
satisfied  for  any  value  of  m.  More  generally,  if  /(c)  is  of  the  form 


fA<)  =  y-,r-gUh) 


(87) 


where  g{x)  is  an  arbitrary  function  of  its  argument,  then  the  compressional  term  can  be 
written 

df 

Cc  —  —Ay-£-.  (88) 

°y 

Thus,  /„(e)  maintains  the  same  form  at  all  times,  with  only  y  varying  in  time  according 
to  dy/dr  —  A2y.  The  parameter  y(t)  is  linearly  proportional  to  the  temperature  of  these 
distributions,  since 


kT  =  -(E) 
<? 


=  y-  [  dxx3/2g(x), 

9  Jo 

and  so  the  temperature  changes  under  compression  at  the  rate 


dT 

dr 


A2(t)T. 


(89) 


(90) 
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If  the  compression  factor  A2  were  independent  of  time,  Cc  would  cause  the 
temperature  of  the  system  to  change  exponentially.  Generally,  however,  A2  depends  on 
time,  and  other  factors  also  contribute  to  /o  in  general.  If  the  compressional  term  is 
combined  with  the  electric  field  heating  terms  (but  not  the  Fokker-Planck  or  inelastic 
terms),  an  equation  for  dyjdr  follows  which  is  very  similar  to  eq.  (81): 

m/2  rm  -  2  -  m(e/y)m/2 


V‘-Yi 


Aim 

z 7^ 


(5) 


3  -  m(e/y)m/2 


(9i : 


As  before,  this  has  self-consistent  solutions  only  when  m  =  5,  when  y(r)  satisfies 

dy  SAj 


dr  2 y3/2 


( 1  +  A2y) 


(92) 


5.  Numerical  method  of  solution 


The  time-dependent  Fokker-Planck  equation,  eq.  (57),  is  numerically  solved,  using 
implicit  finite-differencing.  Implicit  differencing  war  diosen  for  its  stability,  but  the  simpler 
forward-differencing  expression  is  ^iso  used  to  obtain  estimates  of  the  change  in  /  with  a 
particular  time  step. 

For  the  differencing  procedure,  energy  (velocity)  space  is  divided  into  a  possibly  non- 
uniform  grid;  the  grid-point  index  is  denoted  by  a  subscript.  The  energy  at  grid  point  j 
is  t j ,  and  the  grid  spacing  is  given  by  Ae;  =  e;  —  £j_ i.  In  addition,  the  time  evolution  is 
separated  into  discrete  steps  of  duration  At,  with  the  value  of  a  quantity  at  a  given  time 
step  denoted  by  a  superscript.  Eq.  (57)  is  in  the  following  form: 


df  1  d 

Tr=7^Fit)-G[t) 


(93) 


with  F(e)  representing  the  elastic  terms  and  G(e)  representing  the  inelastic  terms.  This  is 
implicitly  finite-differenced  as 


_  j  fn  +  l  -"  ■  ^ 

At 


(/;-  -/;)  = 


rn  +  l/2  _  rn  +  1/2 
j  + 1 12  J-l/2 


v/f7 


Ae 


-  G 


n  +  l/2 


[04] 


j+ 1/2 


with  the  half-gridpoint  values  defined  as  (for  both  space  and  time): 

r*U2  =  \  [/"*'+/"]. 


[9.V 
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Eq.  (94)  is  then  written  out,  with  all  of  the  future-time  quantities  /n+1  on  the  left, 
and  all  the  past-time  quantities  fn  on  the  right-hand  side.  The  abbreviations  M  = 
jV(e)  +  \T^I2dV / dr  and  V  —  D(t)  4-  a£(e )  are  introduced  for  brevity.  The  inelastic 
collision  frequencies  are  denoted  i/*  ■,  where  the  first  subscript  denotes  the  scattering 
process  and  the  second  gives  the  energy  grid-point  index.  Collecting  common  factors, 
the  result  is: 


+  1- 


Ar 


+ 


"Ui-K-i- 
2 


Atj+1  ]+2 


-  —Vn 


Atj  " 


f) 


n+ 1 


4Ae;+^v/i7 


Ar 


At 


+  (‘ + K»  - 


At 


4Aci+iv/€-  [ 


/; 


(96) 


This  can  be  written  in  matrix  form  as 


r^/r1  =  or) 

which  has  the  formal  solution 

fn+1  ={Tn)-lUnfn  (98) 

and  all  of  the  dimculty  in  obtaining  the  solution  lies  in  inverting  the  matrix  T .  If  it  were 
not  for  the  inelastic  collision  sums,  T  would  be  a  tridiagonal  matrix,  whose  inverse  can  be 
rapidly  computed.  With  this  in  mind,  it  is  worthwhile  to  examine  how  important  these 
collision  terms  are.  It  has  been  shown  that  the  overall  contribution  of  inelastic  collision 
to  the  kinetic  equation  is  small.  If  the  distribution  function  does  not  change  rapidly,  the 
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linear  collision  sum  can  be  taken  as  constant  between  time  steps,  and  the  advanced-time 
collision  term  on  the  left-hand  side  of  eq.  (96)  can  be  moved  to  the  right-hand  side.  In 
this  case,  a  tridiagonal  matrix  equation  is  obtained,  which  may  be  written  as  follows: 


cn/j'  +  (i  -  sn)/;+1  +  /T/7-!1  =  Dn  +  /; 


n  rn+1 


(99) 


with 


.4" 

Bn 

Cn 

Dn 


- — -  (Vn  i  ~  —Vn  A 

4Ae,+  iv^7V  ACj  >-h) 


4Ae 


j+iv^ 


Ar 


K»  +  s^») 


4At)+}v^  V 

-cn/;+,  +  bv;  -  4"/;-,  -  Ar^  k/;  -  ^+1(/;+1.) 


(100a) 

(1006) 

(100c) 

(lOOd) 


While  the  implicit  differencing  scheme  is  a  particularly  stable  method  for  the  solution 
of  eq.  (57),  a  simple  forward-differencing  of  the  equation  gives  simpler  results,  from  which 
an  upper  bound  to  the  time  step  At  can  be  derived.  For  the  forward-differencing  of 
eq.  (57),  eq.  (94)  is  replaced  by: 


/n+l  _  rn  _ 

J  J  j  ~ 


At 


Aei+iv/«7 


(/"+.  +  /?)  +  (/"+.  -  /") 


A«,-h 


-  ^  (/;  +  (/;  -  /^,) 


+  (A  r)G] 


(101) 


The  solution  for  the  advanced-time  distribution  function  (corresponding  to  eq.  (99)),  is: 

/;n+1  =  /;  +  2Dn  +  (  At)(7"  (102) 

where  Dn  is  defined  in  eq.  (lOOd). 

This  estimate  of  /"+1  can  be  used  to  fix  an  upper  bound  to  At  such  that  the  change 
over  one  time  step  is  relatively  small.  For  this  estimate,  the  small  inelastic  term  G", 
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involving  sums  over  the  different  present-time  distribution  function  values  /”,  is  ignored. 
If  the  desired  change  is  limited  to  S\ 


K/n+1  -n/r\ <s. 


{ 103) 


then  At  must  satify  the  following  relation: 


.  S  fn 

At  < - 

~  2  (Dn/ At) 


(104) 


where  Dn /At  is  independent  of  At.  In  the  numerical  computations,  (/n+1  —  /")/ /"  was 
limited  to  S  =  0.1  by  means  of  eq.  (104). 

The  procedure  followed  in  the  calculations,  then,  is  to:  (1)  set  up  the  initial  state  by 
specifying  fj  on  the  grid;  (2)  limit  the  time  step  according  to  eq.  (104);  and  (3)  compute 
the  tridiagonal  elements  A ,  B  and  C  of  the  matrix  T  from  the  distribution.  Next,  (4)  the 
matrix  is  inverted  and  (5)  is  used  in  eq.  (98)  to  compute  the  distribution  function  at  the 
next  time  step.  These  steps  are  repeated  to  follow  the  progress  of  the  system  in  time. 

Energy  conservation  and  normalization 

Since  accurate  computation  of  changes  in  the  energy  is  important  in  the  solution,  it  is 
particularly  important  to  properly  difference  the  nonlinear  Fokker- Planck  coefficients  jV(e) 
and  D{e).  It  is  readily  verified  that  the  Fokker-Planck  equation  without  heating  or  cooling 
terms  conserves  energy;  that  is,  the  time  derivative  of  the  kinetic  energy,  as  calculated  by 
eq.  (16),  proves  to  be  zero.  If  this  calculation  is  performed,  it  will  be  seen  that  conservation 
follows  from  the  following  relation: 


dD  r  r 

1T“W7. 


*7(0, 


(105) 


wnich  follows  from  the  definition  of  eq.  (53).  Using  this  formula  rather  than  eq.  (33) 
to  calculate  D  ensures  numerical  energy  conservation.  This  result  was  first  reported  by 
Kho8.  With  this  relation  it  is  simple  to  show  that,  for  a  Maxwellian  distribution  of  the 
form  /(e)  =  Ae~^kT, 


=  kT. 


For  jV(e),  a  straightforward  differencing  of  eq.  (52)  is  sufficient. 


(100) 
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Since  the  energy  of  the  dynamic  system  changes  in  time,  a  normalization  of  the 
energy  coordinate  which  was  based  on  the  system  energy  might  also  change  in  time. 
This  was  certainly  true  of  the  density  normalization,  which  changes  in  time  because  of 
the  compressional  term;  the  density  variation,  however,  could  be  analytically  evaluated 
and  then  incorporated  into  the  compressional  term  (see  eq.  (44)).  This  cannot  be 
done  for  the  energy  variation.  Therefore,  a  time-invariant  normalization  was  employed: 
the  dimensionless  velocity  is  obtained  with  the  initial  thermal  velocity,  not.  the  time- 
varying  one.  At  the  start  of  the  system  evolution,  the  average  normalized  energy  density 
(E)  =  q  =  3/2;  at  later  times  it  varies  freely.  This  means,  for  example,  that  a  later-time 
Maxwellian  has  the  form  /(c)  =  A  exp[-c(kT(0)/kT(t ))).  The  initial  thermal  velocity  is 
simply  a  convenient  scaling  factor,  as  are  the  other  normalization  constants  employed;  of 
course,  they  become  less  cpnvenient  the  more  the  system  changes  from  its  initial  state. 

6.  Numerical  Results 

Relaxation  of  an  initially  Gaussian  electron  distribution  to  a  Maxwell-Boltzmann 
distribution  is  shown  in  fig.  (5).  (Note  that  in  all  plots  of  the  distribution  function,  the 
quantity  plotted  is  the  ratio,  7(c),  of  the  actual  distribution  function  with  a  Maxwellian 
distribution  at  the  same  temperature;  thus,  the  equilibrium  point  in  fig.  (5)  is  a  straight 
line  at  7(c)  =  1.)  I11  this  calculation,  only  the  electron-electron  collision  terms  were 

nonzero.  The  initial  distribution  was  given  by  /(e)  =  a0exp(  — 10[(£  —  p)/u’]2),  where 
£  is  the  dimensionless  velocity  of  eq.  (38)  and  ao,  p  and  w  are  constants  determining 
normalization,  gaussian  center  position  and  width,  respectively.  For  fig.  (5),  the  values 
p  =  0.3  and  w  =  0.3  have  been  taken.  The  result  appears  to  agree  with  that  of  Klio8  for 
the  same  parameters.  It  can  be  seen  in  the  plot  that  the  equilibration  time  at  any  energy 
is  roughly  proportional  to  the  electron-electron  collision  time. 

When  there  is  ohmic  heating,  the  shape  of  the  evolving  distribution  depends  on  (he 
strength  of  the  magnetic  field.  Strong-magnetic-field  ohmic  heating  results  in  a  depressed 
tail,  as  is  shown  in  the  time-development  plots  for  both  weak  heating  (fig.  (6),  m  =  2.5)  and 
strong  heating  (fig.  (7),  m  =  4.0).  These  distributions,  evolving  from  an  inital  Maxwellian. 
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rapidly  take  on  the  self-similar  shape  for  heating  which  was  described  above.  (The  ratio 
of  the  self-similar  distribution,  eq.  (78),  with  a  Maxwellian  has  a  peak  at  an  energy  which 
is  easily  shown  to  be  proportional  to  the  temperature;  the  heating  up  of  the  distribution 
can  be  seen  in  the  figures  by  the  advance  in  the  peak  locaton.) 

Weak-magnetic-field  ohmic  heating  is  centered  at  an  electron  energy  emai  (eq.  (68)) 
which  depends  on  ft-2/3;  the  distribution  is  enhanced  relative  to  the  same-temperature 
Maxwellian  in  the  vicinity  of  €max,  and  high-energy  runaway  production  is  possible.  In 
fig.  (8),  the  time-evolution  of  an  initially  Maxwellian  distribution  is  shown  for  olnnic 
heating  focused  at  e  =  5  (ft  =  0.2324 Z).  As  long  as  the  magnetic  field  is  finite,  the 
highest-energy  electrons  will  always  be  confined,  and  the  distribution  tail  must  be  depleted 
relative  to  a  Maxwellian;  for  low  enough  fields,  however,  an  enhancement  can  be  seen  over 
a  significant  fraction  of  the  high-energy  electron  population. 

Inelastic  collisions  cause  a  relative  enhancement  in  the  distribution  function  tail.  This 
can  be  seen  in  fig.  (9),  which  shows  the  evolution  of  a  system  under  the  sole  influence  of 
inelastic  collisions.  For  simplicity,  only  one  inelastic  process  was  included,  with  a  collision 
strength  ftaj,  =  1;  the  same  behavior  holds  true  in  general.  There  is  a  depletion  in  the 
distribution  function  near  the  threshliold  energy  eai,  but  everywhere  else  the  distribution 
function  is  enhanced  over  the  Maxwellian  value.  This  is  a  surprising  result,  since  the  effect 
of  inelastic  collisions  is  to  move  electrons  from  higher  to  lower  energy  states.  Collisions 
lower  the  temperature  much  more  effectively  than  they  reduce  the  high-energy  electron 
population,  however,  and  the  enhanced  tail  of  fig.  (9)  is  the  result. 

The  effect  of  compression  is  shown  in  fig.  (10),  where  a  non- Maxwell! an  distribution 
(obtained  as  in  fig.  (8)  by  ohmic  heating)  is  compressed.  The  result  is  a  more  or  less 
isomorphic  heating  of  the  distribution.  Just  as  in  fig.  (6)  and  the  following  figures,  the 
peak  in  the  ratio  /o(«)/ /mb  moves  to  the  right  as  the  distribution  is  heated.  As  was  shown 
here,  if  the  initial  distribution  has  the  form  of  the  self-similar  solution  (eq.  (78)),  including 
the  Maxwellian,  compressional  heating  is  completely  isomorphic:  only  the  temperature 
is  changed  by  compression  (with  appropriate  shifts  in  the  peak  of  the  ratio  with  the 
Maxwellian). 

These  figures  illustrate  the  individual  effects  of  the  various  Fokker-Pianck  terms  of  the 
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electron  distribution.  In  a  z-pinch  implosion,  these  terms  are  interacting  and  competing  to 
produce  a  dynamically  changing  distribution  function,  which  may  at  times  be  Maxwellian 
and  at  other  times  be  strongly  non- Maxwellian.  The  numerical  code  spelled  out  here  can 
be  used  both  to  follow  the  changes  of  this  distribution  function  and  to  estimate  the  effects 
of  its  variation.  The  determining  factors  are  the  external  parameters  of  temperature, 
degree  of  ionization,  mean  velocity,  and  electric  and  magnetic  fields.  These  parameters 
could  obtained  from  a  hydrodynamic  model,  for  example,  which  solves  for  the  spatial 
variation  of  a  given  system,  and  the  calculated  velocity  distribution  function  could  be 
used  to  estimate  runaway  production,  radiation  output,  and  the  self-consistency  of  the 
hydrodynamic  model.  Results  from  this  diagnostic  use  of  the  Fokker-Planck  solution  will 
be  reported  in  subsequent  works.  A  hybrid  hydrodynamic-kinetic  model  could  also  be 
constructed,  where  the  transport  coefficients  are  calculated  from  the  true  distribution 
function  instead  of  an  assumed  Maxwellian.  This  procedure  could  be  both  more  accurate 
and  more  robust  in  analyzing  the  variety  of  possible  z-pinch  implosions. 
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APPENDIX  I 

The  Fokker-Planck  Collision  Term 

Here,  the  Fokker-Planck  collision  term  which  is  written  out  in  eq.  (20a)  is  derived. 
This  derivation  is  not  from  “first  principles”;  it  is  mainly  intended  to  enable  comparison 
with  other  forms  of  the  collision  term  found  the  literature.  A  more  detailed  treatment  of 
the  statistics-related  Fokker-Planck  equation  can  readily  be  found  in  the  literature,  along 
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with  the  application  to  Coulomb  scattering  employed  here4. 

We  consider  the  effect  of  collisions  with  one  particle  species  on  the  distribution  function 
of  another  species.  Species  a  has  mass  m  and  distribution  function  f(v),  while  species  b 
has  mass  M  and  distribution  function  F{v).  Changes  are  computed  which  result  from  a 
particles  (e.g.,  electrons),  with  velocity  v,  being  incident  on  b  particles  (e.g.,  ions),  with 
velocity  t»*.  The  relative  velocity  is  t?  =  v  —  v'. 

The  “Landau  form”  of  the  Fokker- Planck  collision  term  is 


YabV.T(v) 


(I.D 


where 

and  Yab  is  defined  in  eq.  (21).  The  gradient  operator  V  acts  on  the  variable  v  unless 
otherwise  specified. 

To  rewrite  eq.  (1.1),  it  is  useful  to  note  that 


t?3 


VvVv# 


(1.3) 


By  commuting  the  V„  operations,  integrating  by  parts  with  V,.>,  and  using  V„.  —  -V,,, 
it  can  easily  be  shown  that 


where 


(tO  =  +  jVV  :  |5(v)/(v)] 


(1.4) 


r/-N  m  +  M  ^  f  ,  F(v  ) 

A(v)  -  ra6V  dv  _ 

M  J  |u  -  v  '| 

(1-5) 

B  =  raiVV  j  dv'\v-  v'\ F(v') 

(Ui) 

The  quantity  A{v)  is  called  the  coefficient  of  dynamic  friction,  and  B(v)  is  called  the 
diffusion  coefficient.  The  integrals  /  du'F(r')d  and  /  dv'F(v')/#  in  the  above  expressions 
are  called  the  Rosenbluth  potentials. 
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The  Cartesian  tensor  expansion  of  eq.  (13)  can  now  be  t\.  cloyed  in  eq.  (1.4)  to  obtain 
equations  which  are  satisfied  by  each  expansion  term.  A  spherical  coordinate  system  is  used 
to  evaluate  the  integrals.  Of  course,  only  a  magnitude- of- v  dependence  will  be  left  in  the 
final  expression,  though  the  direction  of  fi(v)  will  preserve  (to  first  order)  the  directionality 
of  the  system.  It  is  straightforward  to  calculate  the  equation  for  /0(u).  Higher-order 
terms  are  more  complicated,  since  these  involve  vector,  tensor,  etc.  components,  and  the 
coefficients  A(t»)  and  B(v)  must  be  expanded.  This  will  not  be  done  here  (see  Shkarofsky4, 
section  7-5,  for  details). 

Expressions  are  given  below  for  the  required  derivatives  in  spherical  coordinates,  when 
there  is  only  a  |v|  dependence  (with  unit  basis  vectors  {ej}). 

-  t’;  dd> 

V<t>  =  et-^  (1.7) 

v  dv 


Vrr7  j.  -  -  fc  1  d<t>  ,  ViVj  (d24>  ld(t>\ 


=  I<p(v)  -I-  vvxl'(v) 

VV  :  VV<f>  =  V  •  [v  •  VV<£ 

d  ( Vj  {dtp  2  cty  1  1 
dv{  \  v  [  dv  dv  ) 

2  dtp  d2p  2  4  d4'  d24' 

v  dv  dv2  v2  ^  v  dv  dv 2 

In  evaluating  A(v)  and  B(v),  the  following  easily  derived  identities  are  useful: 


d6  sin# 


v—  y'|  max(v,v') 


r  \2(v+ V-)  ’  if  v' <  r; 

/  dB  sin#|v  —  u'|  =  <  ) 

1  M^b)’  «"■>'■ 


(1.10) 


(1.11) 


(1.12) 


(1.13) 


Using  these  Identities  it  is  easy  to  show  that 


/  =47r0/  dv'v'2F0(v')  +  dv'v'F0(v')} 

J  dv'\v  -  v'\Fq{v)  =  J  dv'v'2F0(v' )  +  —  J  dv'v'*Fo(v') 

+  J  dv'v'3 Fo(v')  +  ~ 
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and  so 


An  =  — 


8  ttY. 


V- 


f  dv'v'2 

J  0 


Fo(v') 


(1.16) 


Bo  =  InY"  ( /  [  -  [  dv'v'2  Fo(v')  -  — -j  f  dv'v'*F0(v')  +  \  f  dv'v'F0(v') 
(  [  V  Jo  <51’  Jo  •*  J  v 

+  vv  J  dv'v'*  f(v')  -  -  J  dv'v'2  f(v')  | 


(1.17) 


These  expressions  can  now  be  substituted  into  eq.  (1.4),  which  to  zeroth  order  is 


(t\  = 


-A0{v)f0(v)  +  -V  •  (Bo(v)fo(v)) 


d_ 

dvi 


(Niv)fo(v)  + 

where  the  quantities  N(v)  and  D(v)  are  defined  as 
N(v)=4n  (  dv'v'2Fo(v') 

Jo 


D(v)  =  —  -  /  dv'v'*F0(v')  +  v2  [  dv'v'Fo(v') 
3  v  Jo  Jv 


Finally,  since 


ft*, 
dv 


-  ^ 


r2  dv  ’ 


then  the  lowest-order  expression  for  the  Fokker-Planck  collision  term  is 


( Hl\  _  nY  ° 

V  St  )c  V2  dv 


N(v)Mv)  + 


This  is  just  eq.  (20a). 


(1.18) 

(1-19) 


(1.20) 

(1.21) 

(1.22) 


(1.23) 


APPENDIX  II 

Cartesian-Tensor  Component  Equations 

In  this  section,  component  equations  for  the  first  two  terms  of  the  Cartesian  tensor 
expansion  of  f(r,v,t)  are  derived  from  eq.  (7),  the  Boltzmann  equation  in  stationary 
coordinates.  These  differ  only  slightly  from  eqs.  (18),  which  were  derived  from  eq.  (12). 
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the  Boltzmann  equation  in  moving  coordinates,  and  are  identical  if  V(r)  =  const.  Although 
the  expansion  relies  on  spherical  symmetry,  we  will  use  Cartesian  components  in  t lie 
development.  The  Boltzmann  equation  is 

^+fV/  +  (a  +  uw>VJ  =  C  (II. 1) 

at 

The  expansion  used  is 

f(v)  =  /o  +  -v  •  h 

v 

=  /o  4-  f x  sin  9  cos  <f>  4-  fy  sin  6  sin  <t>  4-  /.  cos  0  (II. 2) 


The  chain  rule  gives  for  V„: 

8v~d_  89_d_  d±d_\  . 

dvi  dv  ^  dv,  89  dv,  d<$>  J  * 

v  8  1  dv  8  1  dv  8 

v  dv  v2  89  89  +  v2  sin2  9  8<t>  8<f> 


( II -3 ) 

(H.4) 


When  the  expansion  of  eq.  (II. 2)  is  substituted  into  the  Boltzmann  equation,  eq.  (II. 1), 
the  various  terms  become: 


df  dfo  , dfx  8fy  .  .  dfz 

m  =  aT  +  ~m  s,n  1 6  cos *  +  W  51,1  * sm  *  +  aT cos " 

v  •  V/  =  v  ( sin  9  cos  <f>,  sin  9  sin  4>,  cos  9)  •  V /0 
4-  v  (sin  9  cos  <f>,  sin  9  sin  <£,  cos  9)  • 


(11.5) 


sin  9  cos  4> V  fx  4-  sin  9  sin  <£V  fy  +  cos  9  V  f. 

. 


(11-6) 


a  •  Vvf  —  ^ 


sin  9  cos  <t> 


dfo  a  jdfx  .  a  .  <  dfy  df. 

— - (-  sin  9  cos  <p— - (-  sm  9  sin<p  — —  4-  cos  9— — 

dv  dv  dv  dv 


4-  -  cos  9  cos  <f>  [fz  cos  9  cos  4>  4-  fy  cos  9  sin  <f  —  /.  sin  9\ 


4 —  sin  <f>  ( fx  sin  <t>  -  fv  cos  <f>] 


1 


ay  |  si 


4-  av  <  sin  9  sin  4> 


4-  sin  9  cos  <t>~~  4-  sin#  sin  4-  cos  6—^ 
dv  dv  dv  dv 


1 


H —  cos  9  sin  <p  \  fz  cos  9  cos  4>  +  fy  cos  9  sin  <f>  —  f-  sin  9] 
v 
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4-  ^  cos  4>  [~fx  sin  4>  4-  fy  cos  4>\  j 

4-  az  /  cos  9  4-  sin  9  c  os<p~~-  4-  sin  9  sin  d>~l—  4-  cosB-J- 

"  )  dv  dv  0  v  dv 

4-  -  sin#  \  —  fx  cos  9  cos  4>  -  /„  cos  #  sin  <£  H-  /-  sin  0]  > 

v  J 


v  x.  uj  •  V,,/  =  —  u>  •  y  x  Vvf 


=  ljx  [—fy  cos  9  4-  fz  sin  9  sin  4>\  4-  u>y  [ fx  cos  9  -  fz  sin  9  cos  <$>\ 
4-  u>2  [sin  9(  -fx  sin  <t>  +  fy  cos  <£)] 


(II.8) 


To  obtain  an  equation  for  /0,  the  Boltzmann  equation  resuting  from  these  substitutions 
is  integrated  over  solid  angle  (d2Q  =  sin  9  d9d<j>).  The  results  are 


/f?S  in  0^  =  4*^ 

J  v  •  V /  sin  8  d8d<t>  =  1»  V  ■  /i 

/_  -  ,  .  .  ,  4tt  d  ,  _  - .  8tt  _  - 

a  •  V vf  sm 9 d9d<j>  =  —  —(a  •  fi )  4-  —a  •  /j 

3  ar  3v 


( II .9 ) 


(11.10) 


(11.11) 


To  obtain  an  equation  for  /j ,  the  same  procedure  is  followed,  multiplying  the  equation 
in  turn  by  sin<£,  cos<f,  and  cost?.  After  numerical  constants  are  divided  out,  the  resulting 


two  equations  are4 


d/o  /*,  1  0  2-  T  \-n 

~dt+3v'fi  +  z^d:  r  a'h\  ~Cfp  +  c 

^4.-V/,4a^4Jx/:  =  <?, 


1  d 


(11.12) 


(11.13) 


APPENDIX  III 

Electric  field  -  Laser  Heating  equivalence  to  Ohmic  Heating 

While  this  study  is  primarily  concerned  with  ohmic  heating,  there  is  a  direct  analogy 
with  laser  heating  by  inverse  bremsstraiilung,  which  has  been  well  studied  in  the  literature.6 
In  this  appendix,  the  electric  field  heating  term  from  laser  heating  is  derived.  It  is 
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remarkable  that,  under  certain  conditions,  this  is  identical  in  form  to  that  for  ohmic  heating 
derived  earlier. 


The  starting  point  for  the  derivation  is  eq.  (33)  for  the  electric  field  heating  term  in 
the  Fokker-Planck  equation,  and  eq.  (34)  for  f\ .  With  laser  heating,  time  variations  in  f\ 
cannot  be  ignored;  they  are,  however,  periodic.  We  assume  that  any  background  magnetic 
field  is  negligible  compared  to  the  electric  field  of  the  laser  light,  so  u>  =  0  in  eq.  (19b). 
Particles  are  accelerated  by  the  laser  electric  field,  which  has  intensity  I  —  cE2/4n, 
wavelength  A,  and  oscillation  frequency  u >i  =  2irc/X.  The  maximum  oscillation  velocity  of 
the  particles  is  u0,  which  is  related  to  the  intensity  and  frequency  by 


•>  2 

j  _  rn-uzv'c 
4ne2 


(HI.l) 


If  I  is  in  W/ctir,  A  is  in  microns,  and  kT  is  in  keV,  then  the  ratio  of  the  oscillation  velocity 
to  the  thermal  velocity  (with  q  —  3/2)  is  \mvl/kT  =  3.73595  x  10~"16/A2 /kT.  The  laser 
acceleration  is 

a(t)  =  v0u>l  cos ijj it.  (III. 2) 

Under  these  assumptions,  eq.  (34)  for  fi  becomes 

”-  +  t'e,/i  =  -a(0^  (HI*3) 

As  with  ohmic  heating,  it  is  assumed  that  fi  varies  in  time  much  more  rapidly  than 
/o •  But  instead  of  looking  for  equilibrium  solutions  for  fi(v,t),  eq.  (III. 3)  is  solved  for 
the  time-dependent  /lt  and  this  solution  is  averaged  over  a  full  period  of  the  (rapid)  laser 
oscillation.  The  solution  to  the  homogeneous  version  of  eq.  (III. 3)  (i.e.,  a(t)  =  0)  is 


hit)  =  Ae~Vr'i  ( III  .4 ) 

with  A  an  arbitrary  constant.  The  solution  when  a(t)  ^  0  can  be  expressed  in  the  same 
form: 

fi(i)  =  B(t)e~l'-t  (111.5) 

where 

B(t)  =  -  /  dt'.  (III. 6) 

Jtn  9v 
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The  lower  limit  of  integration  is  determined  by  the  boundary  condition,  which  in  this  case 
is  set  by  the  requirement  that  lim  fi(t)  =  0.  To  realize  this,  the  acceleration  a(i )  from 

t  — *  —  OO 

eq.  (III.2)  is  “adiabatically  turned  on”  by  multiplying  it  by  a  factor  exp(d),  where  f  is 
infinitesimal,  for  t  <  0.  When  this  a(t)  is  used,  and  if  dfo/dv  is  assumed  to  be  independent 
of  time  (at  least  on  the  time  scale  of  the  f\  variations),  the  general  solution  for  /i ( / )  is 


/i(0  =  -iTo  — 


ul 


dfo 


* L  +  vei  dl' 


(we,  cos uJit  +  (jJL  sina/£,f) . 


(in.: 


For  the  calculation  of  /0,  the  rapidly  oscillating  values  of  f\  are  averaged.  This  is 
accomplished  as  follows: 

fUJl  t=27T 


2rr  j„lt=0 


Vtivl/2  dfo 


1  +  (i/e,/u>L)2  dv 


( III.8 ) 


If  this  expression  is  used  for  a  ■  /j,  the  laser  heating  term  in  eq.  (33)  is 


ca  =  cL  = 


njY'jvl  d 
6v2  dv 


1 


1 


_ dfo 

v  1  +  (V'Jul)2  dv 


( III .9 ) 


When  expressed  in  terms  of  the  dimensionless-energy  representation  of  eq.  (57)  (fi^  = 
iml/vr),  the  connection  with  the  ohmic  heating  term  can  be  clearly  seen: 


where 


r  «  #<(e) 
Cl  =  77T' 


r  Vo  2?r Z  rec 

oil  =  &  T"  — - “z: l 

v: 


(III. 10) 


th 


q  u>‘LkT 

?  ^3/72^3 


3Z\^1  +nle3/Z2q 


(III. 11a) 
(III.llM 


The  results  for  ohmic  heating  and  for  laser  heating  are  thus  equivalent,  with  an  effective 
laser-heating  “magnetic  field” 

f =  Qeb  (III. 12) 
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and  an  effective  laser  “electric  field” 


p  -  L  v° 

L~  V2  vtk 

where  nmi'2h/2  is  the  thermal  energy  density  in  the  plasma. 

APPENDIX  IV 

Numbers  and  Units 

From  the  1986-87  CRC'  Handbook: 
c  =  2.99792458  x  1010  cm/sec 
m  =  9.109534  x  10~28  g 

e  =  1.6021892  x  10”19  C  =  4.8032423  x  10"10  statcoul 
ft  =  1.0545887  x  10 ~ 2 ‘  erg-sec 
Ry  =  13.605826  eV 
rf  =  2.8179380  x  10-13  cm 

mec2  =  5.1100258x10s  eV  =  8.1872281  x  10"7  erg 
With  these  values,  and  with  kT  measured  in  eV  and  ne  in  cm'3,  the  following  numeric 


formulas  can  be  given: 

(kT)3?2  ' 

log  A  =  23.463306  +  log  - — -L-  (IV.l) 

L  v** 

uR  =  2.1030335  x  10"6  ~r-jr  log  A  sec'1  (IV.2) 

(kT)3'2 

up  =  5.641458  xlO4  ^  (IV. 3) 

Eq  =  9.7l48016x  103  (e/c)upV3kT  statvolt/cm 

=  8.0877853  x  10“ 16  vr  \fkT  statvolt/cm  (IV.4) 

Bd  =  5.685680  x  10-8  vr  Gauss  (IV.5) 
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Figure  4.  Exponent  m  in  the  self-similar  solution  for  the  distribution  function  under 
high-magnetic-field  heating,  as  a  function  of  the  electric  field  in  units  of  ED,  the  Dreicer 
field- 
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Figure  5.  Time-evolution  of  the  distribution  function  from  an  initial  Gaussian  shape 
/(e)  =  a0e~,0!(*-0'3)/,0,3!J  to  the  equilibrium  Maxwellian.  Shown  is  the  ratio  f(e)/fntB(t) 
of  the  actual  distribution  with  a  Maxwellian  at  the  same  temperature.  The  various  curves 
correspond  to  the  times:  A)  initial  condition;  B)  0.02 ;  C)  0.04;  D)  0.08;  E)  0.16;  F)  0.32; 
and  C)  0.64 ,  in  units  of  the  thermal  electron  collision  time  tr. 


Electron  Distribution,  m=2.5  (MB) 


n  1.0  Hk- 


E  * 


*\ - »' 

C* 


>  >  /* 

v  '  \  v  “  < 


X  * 


2  3  4  5  6  7 

Energy  (thermal  units) 


9  10 


Figure  6.  Time-evolution  of  the  distribution  function  from  an  initial  Maxwellian,  under 
the  influence  of  strong-magnetic-field  ohmic  heating.  The  electric  field  is  such  that  the 
exponent  m  =  2.5.  Shown  is  the  ratio  /(€)//a/b(c)  of  the  actual  distribution  with  a 
Maxwellian  at  the  same  temperature,  at  times  A)  initial  condition;  B)  0.50;  C)  1.00;  D) 
2.01;  E)  3.01;  F)  4.02;  G)  5.01;  and  If)  10.06,  in  units  of  the  thermal  electron  collision  time 
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Electron  Distribution,  m=4 . 0  (MB) 
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Figure  7.  Time-evolution  of  the  distribution  function  from  an  initial  Maxwellian,  under 
the  influence  of  strong-magnetic-field  ohmic  heating.  The  electric  field  is  such  that,  the 
exponent  m  =  4.0.  Shown  is  the  ratio  f(e)/fMB{t)  of  the  actual  distribution  with  a 
Maxwellian  at  the  same  temperature,  at  times  A)  initial  condition;  B)  0.06;  C)  0.11;  D) 
0.21;  E)  0.31;  F)  0.40;  G)  0.50;  and  H)  1.00,  in  units  of  the  thermal  electron  collision  time 
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Figure  8.  Time-evolution  of  the  distribution  function  from  an  initial  Maxwellian,  under 
the  influence  of  weak-magnetic-field  ohmic  heating.  Heating  is  focused  at  c  =  5.  The 
electric  field  is  such  that  the  exponent  m  —  2.5.  Shown  is  the  ratio  /(«)//mb(()  of 
the  actual  distribution  with  a  Maxwellian  at  the  same  temperature ,  at  times  A)  initial 
condition;  B]  0.53;  C)  1.08;  D)  2.12;  E)  3.15;  F)  4.09;  and  G)  4.89,  in  units  of  the  thermal 
electron  collision  time  tr. 
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Figure  9.  Time-evolution  of  the  distribution  function  from  an  initial  Maxwellian  under 
the  influence  of  inelastic  collisions  only.  Only  one  inelastic  processes  is  present,  with  a 
threshold  at  e  =  2.  The  collision  strength  has  been  set  to  17  =  1000,  an  unrealistically  high 
value,  for  display  purposes.  Shown  is  the  ratio  /(«)/ /mb(c)  °f  ^,e  actual  distribution  with 
a  Maxwellian  at  the  same  temperature,  at  times  A)  initial  condition;  B',  0.10;  C)  0.29;  D) 
0.48;  and  E)  0.67,  in  units  of  the  thermal  electron  collision  time  tr. 
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Figure  10.  Time-evolution  of  a  distribution  function  from  an  olimically- heated 
distribution  under  the  influence  of  a  constant  compression ai  heating  only.  (Initial 
distribution  is  that  of  fig.  (8)  at  r  =  1.078).  For  this  plot,  the  compression  coefficient 
was  set  to  f  t^V  •  V  =  1  (see  eq.  (57)).  The  different  curves  correspond  to  times  A)  initial 
condition;  B)  0.11;  C)  0.21;  D)  0.31;  E)  0.41;  and  F)  0.92,  in  units  of  the  thermal  electron 
collision  time  tr. 
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